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We report on the realization of a quantum degenerate atomic Fermi gas in an 
optical lattice. Fermi surfaces of noninteracting fermions are studied in a three- 
dimensional lattice. Using a Feshbach resonance, we observe a coupling of the 
Bloch bands in the strongly interacting regime. 



1. Introduction 

The exploration of quantum degenerate gases of fermionic atoms is driven 
by the ambition to gain deeper insight into long-standing problems of quan- 
tum many-body physics. So far, however, the analogy to an electron gas 
in a solid is limited since there the electrons experience a periodic lattice 
potential. The lattice structure is in fact a key ingredient for most models 
describing quantum many-body phenomena in materials. We access this 
regime by preparing a degenerate atomic Fermi gas in the crystal structure 
of an optical lattice. 

In our experiment we load a noninteracting gas of '**'K atoms into a 
three-dimensional optical lattice with simple cubic symmetry and directly 
image the Fermi surfaces. Gradual filling of the lattice transforms the 
system from a normal state into a band insulator. Previous experiments 
with far-detuned three-dimensional optical lattices were always carried out 
with bosonic atoms^'^'^, and experiments with fermions were restricted to 
a single standing wave^. 

Interactions in the atomic Fermi gas can be tuned by using a Feshbach 
resonance between different spin components of the gas. By ramping the 
magnetic field into vicinity of the resonance, we increase the interactions 
between two particles residing on the same site of the lattice. In this manner 
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we dynamically induce a coupling between the lowest energy bands. 

2. Preparing the degenerate Fermi gas in the 3D optical 
lattice 

The apparatus and the procedure we use to create a degenerate atomic 
Fermi gas is described in previous work-'^. In brief, we use bosonic *^Rb to 
sympathetically cool a spin-polarized gas of fermionic '^''K atoms to tem- 
peratures of T « 0.3 Tp (Tp = 260 nK is the Fermi temperature of the 
noninteracting gas). The potassium atoms are then transferred from the 
magnetic trap into a crosscd-beam optical dipole trap with a wavelength 
of A = 826 nm. There we prepare a spin mixture with (50 ±4)% in each of 
the |F = 9/2, niF = -9/2 > and |F = 9/2, uif = -7/2 > spin states using 
a sequence of two radio frequency pulses. After further evaporative cooling 
in the optical trap, we reach temperatures between T = 0.2 Tp and 0.25 Tp 
with 5 X lO** to 2 X 10^ particles, respectively. 

Prior to loading the atoms into the optical lattice we tune the magnetic 
field to B = (210 ± 0.1) G, such that the s-wave scattering length between 
the two states vanishes. We exploit the magnetic Feshbach resonance be- 
tween the |F = 9/2, mp = -9/2 > and |F = 9/2, tof = -7/2 > states'^ 
which is centered at Bq = 202.1 G and has a width of AS = 7.8 G. In this 
way we produce a two-component Fermi gas without interactions. 

Then the standing wave laser field along the vertical z-axis is turned 
on. Subsequently, the optical dipole trap along the y-axis is turned off and 
a standing wave laser field along the same axis is turned on, followed by 
the same procedure along the a;-axis. In order to keep the loading of the 
atoms into the lattice as adiabatic as possible, the intensities of the lasers 
are slowly increased (decreased) using exponential ramps with durations 
of 20 ms (50 ms) and time constants of 10 ms (25 ms), respectively. In its 
final configuration the optical lattice is formed by three orthogonal standing 
waves with a wavelength of A = 826 nm, mutually orthogonal polarizations 
and 1 / e^-radii of 50 /um (cc-axis) and 70 fiuv (y-axis and z-axis) . The lattice 
depth Vo is calibrated by modulating the laser intensity and studying the 
parametric heating. The calibration error is estimated to be < 10%. 

3. Observing the Fermi surface 

The potential created by the optical lattice results in a simple cubic crystal 
structure. The Gaussian intensity profiles of the lattice beams give rise to 
an additional confining potential which varies with the laser intensity. As 
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a result, the sharp edges characterizing the T = distribution function for 
the quasimomentum in the homogeneous case are expected to be rounded 
off. Noninteracting fermions on a lattice with an additional harmonic con- 
finement are described by the tight-binding Hamiltonian^ 

H = -JY1 ^J^^' + ^ XI i^ext,xXi + '^Lt,yyi + ^ext.z^f ) ^« • (1) 
{hJ) » 

Here cj and Cj denote the fermionic creation and annihilation operators for 
a particle at site ?', and n; = cjc,; is the mrniber occupation of site i. The 
first term with the tunneling matrix element J corresponds to the kinetic 
energy the particles gain through delocalization over the neighboring lattice 
sites. The second term accounts for a site-specific energy offset due to 
the harmonic confinement, which is represented by its trapping frequencies 
'^ext.a (with a = x,y,z). As the values of Wext depend on the intensity of the 
lattice, changing the tunnel coupling J by varying the depth of the periodic 
potential results in a simultaneous change of the external confinement. 

The harmonic confinement leads to a curvature of the Bloch bands in 
position space and, as a consequence, spatial regions with different fillings. 
For a given Fermi energy , the filling in the center is highest and falls off 
to zero to the outer regions of the trap. When Ep exceeds the band width, a 
band insulating region appears first in the center of the trap, while regions 
farther out still have a filling below unity. However, due to the curved 
bands in position space, these atoms may still be spatially localized to a 
few lattice sites for small band widths'^. 

In the experiment we probe the population within the Brillouin zones 
by ramping down the optical lattice slowly enough for the atoms to stay 
adiabatically in the lowest band whilst quasimomentum is approximately 
conserved^. We lower the lattice potential to zero over a timescale of 1 ms. 
After an additional 1 ms we abruptly switch off the homogeneous magnetic 
field and allow for a total of 9 ms of ballistic expansion before we take an 
absorption image of the expanded atom cloud. The momentum distribution 
obtained from these time-of-flight images reproduces the quasimomentum 
distributions of the atoms inside the lattice. 

With increasing particle number and confinement, the filling of the lat- 
tice sites increases. The initially circular shape of the Fermi surface (fig- 
ure la) develops extensions pointing towards the Bragg planes (figure lb) 
and finally transforms into a square shape completely filling the first Bril- 
louin zone deeply in the band insulator. We have observed population of 
higher bands if more atoms are filled into the lattice initially (figure Ic). 
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Figure 1. Development of the Fermi surface in quasimomentum space when more and 
more fermions are filled into the lattice, a) 3500 atoms per spin state and Vb = 5 Er- 
b) 15000 atoms per spin state and Vb = 5Er- c) 100000 atoms per spin state and 
Vb = 20 Et- When the Fermi energy exceeds the initial width of the lowest Bloch band, 
higher Brillouin zones begin to become populated. As the band structure of the second 
and the third band in a 3D optical lattice overlap partially, both the second and the third 
Brillouin zone are being filled with atoms. Each of the three displayed quasimomentum 
distributions is obtained by averaging and smoothing 4 absorption images taken with 
the same parameters. 



4. Tuning the interactions 

By using a Feshbach resonance, it is possible to tunc continuously from 
attractive to repulsive s-wave interactions. When the on-site interac- 
tion is much smaller than the band gap, the physics in the lowest Bloch 
band can be described in the tight-binding approximation by a Hubbard 
Hamiltonian-'^'' . The phase diagram of the Fermi-Hubbard Hamiltonian is 
very rich, including Mott insulating^^ and antiferromagnetic phases^^. We 
explore the limit where the wells of the lattice can be regarded as an array 
of isolated harmonic oscillator potentials, and we use a Feshbach resonance 
to tune the interactions between the particles in each of the wells. 

In deep optical lattice potentials where tunneling is negligible, the Bloch 
bands are nearly flat and can be approximated by the levels of the harmonic 
oscillator inside each lattice well. When the s-wavc interaction is changed 
on a time scale short compared to the tunneling time between adjacent 
potential minima, we may regard the band insulator as an array of more 
than ten thousand independent harmonic potential wells, each of which 
is populated by two particles. Solving the Schrodinger equation of two 
particles with a contact interaction in a harmonic oscillator yields^^ 

^ / /^ r(-i;/(2M + 3/4) \-^ 

flho V n-E/{2hoj) + l/4)J ^ ' 

where E denotes the eigenenergies for a given scattering length a^, and u) 
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is the trap frequency. The osciUator length aho = yJh/{mLo) for a single 
lattice well in our system is on the order of 1300 oq. Figure 2a displays the 
energy spectrum as a function of the scattering length a^. 




Scattering length [a^o] Magnetic field [G] 



Figure 2. a) Energy spectrum of two strongly interacting particles in a three- 
dimensional harmonic oscillator potential which is given by equation (2) in the center-of- 
mass frame. The asymptotic energies at as — » itoo are indicated by the dashed lines. The 
arrows show the direction of our sweep over the Feshbach resonance, b) Experimental 
method to observe interaction-induced transitions between Bloch bands. Two Feshbach 
resonances between the \F = 9/2, mp = —9/2 > and \F = 9/2, mp = —7/2 > states 
(dashed line, left) and the \F = 9/2,mp = -9/2 > and \F = 9/2, mp = -5/2 > states 
(dotted line, right) are exploited to tune the interactions in the gas. Using an rf pulse, 
we transfer the atoms from the \F = 9/2, nip = -7/2 > to the \F = 9/2, mp = -5/2 > 
state prior to sweeping over the Feshbach resonance at 224.2 G. 



On the Feshbach resonance, the scattering length has a pole {as — > ±oo). 
When sweeping adiabatically over the Feshbach resonance from magnetic 
fields below the resonance, we first increase the repulsive interaction be- 
tween the two atoms in the oscillator ground state from the value of the 
background scattering length to —t +oo on the resonance. There the en- 
ergy level of the lowest harmonic oscillator level becomes degenerate with 
the energy of the second next level, which is lowered from its noninteracting 
value because of the divergence —t — oo on the resonance. If projected 
onto the noninteracting eigenstates, this corresponds to a superposition of 
many higher oscillator states. When continuing the magnetic field sweep 
to values where the scattering length is again near zero, the atoms may 
adiabatically follow the higher oscillator level. This simple picture of two 
strongly interacting atoms in a well gives a qualitative understanding on 
what happens when atoms in a three-dimensional optical lattice cross a Fes- 
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hbach resonance: The higher energy bands get populated due to interaction 
induced coupling between the Bloch bands. 

We experimentally investigate the interacting regime in the lattice start- 
ing from a noninteracting gas deep in a band insulator with Vx = 12 
and Vy = Vz = 18 E,- and corresponding trapping frequencies of lu.,. = 
27r X 50 kHz and Uy = Uz = 2tt x 62 kHz in the individual minima. A 
short radio-frequency pulse of 40 fis is applied to transfer all atoms from 
the |F = 9/2, mp = -7/2 > into the \F = 9/2, mp = -5/2 > state, with 
the atoms in the \F — 9/2, mp = —9/2 > remaining unaffected. Then 
we ramp the magnetic field to different final values around the Feshbach 
resonance^^ located at Bq = 224.2 G which has a width of AB = 9.7 G. 
The sweep approaches the Feshbach resonance from the side of repulsive 
interactions and crosses the resonant point where the scattering length di- 
verges towards the side of attractive interactions, as depicted in figure 2b. 
When using this direction of the sweep there is no adiabatic conversion to 
molecules. After turning off the optical lattice adiabatically and switching 
off the magnetic field we measure the momentum distribution. To see the 
effect of the interactions we determine the fraction of atoms transferred into 
higher bands. 

Figure 3a shows the measurement results when ramping the magnetic 
field to different final values around the Feshbach resonance with a fixed 
inverse sweep rate of 12/is/G. For final magnetic field values well above 
the Feshbach resonance we observe a significant increase in the number 
of atoms in higher bands along the weak axis of the lattice, demonstrat- 
ing an interaction induced coupling between the lowest bands. When we 
choose symmetrical lattice intensities, the atoms in the higher bands are 
distributed over all three dimensions since the oscillator levels are degener- 
ate. By using a single, weaker axis we lift this degeneracy and transfer only 
to the higher bands along this direction, which increases our signal-to-noise 
ratio crucially. 

Initially, one might expect that all atoms should be transferred to the 
higher bands. However, the process could be limited by the actual number 
of doubly occupied lattice sites which constitutes an upper bound for the 
transferred fraction. The inhomogencous filling of the lattice and the finite 
temperature give rise to a significant amount of singly occupied sites, espe- 
cially in the regions away from the center. When we vary the hold time at 
the final magnetic field far above the resonance, we do not observe a sig- 
nificant dependence on the measured fraction of atoms in the higher bands 
(figure 3b). Very recent theoretical studies^^'^^ are beginning to shed light 
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Figure 3. Interaction induced transition between Bloch bands, a) A fraction of the 
fermions is transferred into higher bands using a sweep across the Feshbach resonance 
(filled symbols). The line shows a sigmoidal fit to the data. The open symbols show a 
repetition of the experiment with the atoms prepared in the spin states \F = 9/2, mp = 
—9/2 > and |_F = 9/2, mp = —7/2 > where the scattering length is not sensitive to 
the magnetic field. The magnetic field is calibrated by rf spectroscopy between Zeeman 
levels. Due to the rapid ramp the field lags behind its asymptotic value and the horizontal 
error bars represent this deviation. The data is taken from reference [1] . The inset shows 
the quasimomentum distribution for a final magnetic field of 233 G, for which the images 
of 4 measurements have been averaged and smoothed. Arrows indicate the atoms in 
the higher bands, b) Measured fraction in the higher bands versus the hold time at 
233.8 G after the sweep. The vertical error bars show the statistical error of 4 repetitive 
measurements. 



onto physics beyond the single-band description of the Feshbach resonance 
in the lattice. 



5. Conclusions 

In our experiments, we have realized a degenerate Fermi gas with tunable 
interactions in a three-dimensional optical lattice. Wc have demonstrated 
the control over parameters of the system such as filling and interactions. In 
the noninteracting, static regime, we image the shape of the Fermi surface 
for different characteristic densities. Our measurements in the strongly in- 
teracting regime where we observe a dynamical coupling of the Bloch bands 
pose challenges for the present theoretical understanding of many-particle 
fermionic systems in optical lattices. This very generic implementation of a 
fcrmionic many-particle quantum system on a lattice is expected to provide 
new avenues to intriguing phenomena of quantum many-body physics. The 
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unique control over all relevant parameters in this intrinsicially pure system 
allows us to carry out experiments which are not feasible with solid-state 
systems. 
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